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Abstract 

We review some recent results on recursion relations which help evaluating arbitrary non-diagonal, radial hydrogenic 
matrix elements of T" and of pr^ ((3 a Dirac matrix) derived in the context of Dirac relativistic quantum mechanics. 
Similar recursion relations were derived some years ago by Blancliard in the non relativistic limit. Our approach is 
based on a generalization of the second hypervirial method previously employed in the non-relativistic Schrodingcr 
case. An extension of the relations to the case of two potentials in the so-called unshifted case, but using an arbitrary 
radial function instead of a power one, is also given. Several important results are obtained as special cases of our 
recurrence relations, such as a generalization to the relativistic case of the Pasternack-Sternheimer rule. Our results 
are useful in any atomic or molecular calculation which take into account relativistic corrections. 

1 Introduction 

In atomic or molecular physics it is customary the use of sum rules to handle the matrix 
elements that appear in many computations pQ - [Oj. One of the conspicuous examples is the 
Blanchard relation ^ which is a useful recurrence formula for non diagonal, non-relativistic 
arbitrary matrix elements, of the form (ni/i|r^|n2/2)- Here \nl) stand for non-relativistic 
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hydrogenic radial energy eigenstates, and A is an arbitrary, even non- integer- power of r} 
According to this relation, once we know any three successive matrix elements of powers 
of the radial coordinate, r, any other of these elements can be deduced in terms of the 
three previous ones. The Blanchard recurrence relation was derived more thirty years ago 
using a calculation-intensive method. In the intervening years a simplified method has been 
proposed for deriving the Blanchard rule and, also, new potentially useful relations [H] - |12j . 
Different approaches have been studied for obtaining sum rules between hydrogenic matrix 
elements of these sort. Some of them are based on the non relativistic hypervirial theorem 



Most of the sum rules yet discovered are of a non relativistic nature despite the physical 
and chemical interest for obtainig relativistic results P, (TUI, [23, [2^1 • Nevertheless, there 
has been some effort in constructing new sum rules in the relativistic and quasi-relativistic 
approach jT2]- ^1], [21], [23 • In this review we want to present some of the results obtained 
by us in the relativistic case, using an approach inspired on the non relativistic hypervirial 
method proposed by some of us some years ago [E], [ZI- 

In atomic and molecular physics the non-relativistic approach usually means non rela- 
tivistic matrix elements of powers of a radial coordinate between states of the system at hand 
[H] - [21] ■ But in the relativistic approach we must to take into account, in the general case, 
the presence of the Dirac's (3 matrix in an explicit way in our calculations. This point seems 
natural since the /3 matrix has an eigenvalue -|-1 for the positive energy eigenstates of the 
Hamiltonian and —1 for the negative ones, giving some diferences in the treatment of both 
matrix elements. This point gives in general, two series of sum rules, one for matrix elements 
without the (3 matrix and other one for matrix elements with the f3 explicitly appearing. To 
this end we employ a relativistic calculation inspired directly on the hypervirial method [H] 
to deduce new recurrence relations for the, in general, non-diagonal radial matrix elements 
of succesive powers of and of jSr^ — where /5 is a 4x4 Dirac matrix JH] — for relativistic 
hydrogenic states in the energy basis. The assumptions we use here are that the nucleus is 
point-like and fixed in space, and that a description using the Dirac equation is, of course, 
valid. We first study the recurrence relations in the general case, in which the matrix ele- 
ments are taken between states with different principal quantum numbers rii ^ n2,, different 
total angular momentum quantum numbers ji ^ j2, rrij^ ^ m^j, and different parity. For 
the sake of convenience, we employ the quantum number e = (— instead of parity 
for labelling the hydrogenic eigenstates, where ei 7^ €2 is defined by 



As we mentioned above, we find that in general the recurrence relations depend on 
matrix elements of both powers of r and of f3r, in practical terms this means that we need 
two recurrence relations as the relativistic version of the single-equation Blanchard relation. 
Given its special interest, we in particular study the case where the total angular momentum 
and parity become equal, ji = j2 and ei = €2, in the two states — not mattering the relative 

^Though perhaps of relatively little physical interest, it is nevertheless worth noting that the exponent, 
A, can be complex, see [H] for details 



[SI, [3, cni, m- 
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values of the principal quantum number n. We also address the completely diagonal case 
where ni = n2, ji = j2, and ei = €2- Both of the particular cases mentioned above require 
special treatment for avoiding possible divisions by zero in the general expressions; such 
results are immediately used to obtain a relativistic version of the Pasternack-Sternheimer 
relation and to obtain an expression for the relativistic virial theorem ^2); EI]- 

The results reported in this paper are important because the link between quantum cal- 
culations and experimental results is made at the level of matrix elements of the appropriate 
operators. In atomic and molecular physics this usually means matrix elements of powers 
of a radial coordinate between states of the system at hand [20] - - [2H|; 112] • But 

matrix elements of more general radial functions are also very useful |2S] - EE]- Any con- 
tribution making less cumbersome the evaluation of a series of these elements is potentially 
very useful. In nonrelativistic quantum mechanics the importance of hypervirial results and 
other related techniques follows from this fact, since the task of calculating matrix elements 
can be indeed simplified These techniques are also important for atomic physics 

in the relativistic realm J2]- CS]- This can be significative at present given the precision 
attained in atomic physics experiments using synchroton radiation which is typically less 
than the expected values of relativistic corrections in some processes jHT] - |¥T] . 

In section 2 of this paper, we review the non-relativistic second hypervirial scheme and 
we use it to derive the Blanchard relation. In section 3, we deduce a radial Hamiltonian, 
completely equivalent to the Dirac equation for central potentials and we employ our result 
to implement a corresponding hypervirial result in relativistic quantum mechanics, and we 
proceed to use it to deduce by a long, but direct calculation, the relativistic recurrence 
formulae. In section 4 we study in particular the diagonal case (j^ = ji, ^2 = ei), which 
needs a special treatment to avoid division by zero, to derive the relativistic Pasternack- 
Sternheimer rule and use it (when ni = ^2) to obtain a version of the relativistic virial 
theorem. In section 5 we use Dirac equation to obtain analytic expressions for the bound- 
bound, diagonal and non-diagonal matrix elements for arbitrary powers of r. As it becomes 
evident, such results are rather cumbersome for relatively large values of the power; for small 
values, on the other hand, they are better regarded as starting values for the recurrence 
relations derived in section 3 and 4 of this article. In section 6 we explore the extension of 
our method to the calculation of generalized recurrence relations for the case of two different 
potentials. We call the results generalized recurrence relations as they relate matrix elements 
of a radial function with those of its derivatives. Also because in the case the function is 
chosen as a power of the radial coordinate and the potential is of the Coulomb type we 
obtain the sum rules of the preceding sections. 

2 The non-relativistic hypervirial method 

In this section we develop an alternative approach based on a generalization of the virial 
method to obtain the original Blanchard relation. We remind that both Blanchard relation 
and its predecesor the Kramers selection rule, were originally obtained employing directly 
the Schrodinger equation together with appropriate boundary conditions in a rather long 
way IE] 5 CP ■ The method proposed originally by us is a much simpler one . It is based on a 
generalized hypervirial result and certain Hamiltonian identities that has been developed to 
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simplify the computations. This technique seemed also to us an appropriate starting point 
for deriving relativistic recurrence formulae, and it is with that point of view that we review 
in this section the hypervirial method as it is applied in non-relativistic quantum mechanics. 
We employ atomic units h = m = e = 1. 

The idea is to start with the radial Schrodinger equation for a central potential V{r) 
written in the form 

Hk\nklk) = En^^ijukh), (2) 

where \nklk) = i^n^iki''^) ^n^h ^'^^ respectively, an energy eigenfunction and its energy 
eigenvalue corresponding to principal and angular momentum quantum numbers, Uk and 
Ik'i k is a. label that we will employ to identify the left bra and right ket in more complex 
epressions. Hk, the non-relativistic radial Hamiltonian, is given by 

H,._l|l-iiL + k(|±i)+VM. (3) 

/ dr'^ r ar zr"^ 

2.1 A nonrelativistic hypervirial relation 

Although we want to calculate the radial matrix elements for terms of the form r^, it is best 
for our purposes to consider first matrix elements of an arbitrary radial function /(r); with 
such a choice we can readily show 

{E, - Ek){n,k \f{r)\nklk) = (r^J. |(-^/" " f'^ - -/' + ^^nkh), (4) 

\ 2 ar r z r'^' 

where we use A^^ = U {U + 1) — (4 + 1)^ Ek = En^i,., and the primes stand for radial 
derivatives. Please recall that the matrix element of an arbitrary radial function /(r) is 

{nili\f{r)\nklk) = r'^'ilj*r{r)f{r)^n^i^{r)dr. (5) 







To establish the result we are after, we apply the previous result to the radial function 
e(r) = HJir) - fir)Hk, to find 

2{Ei-Eky{nih\f{r)\nklk) = 
{n,k\{H,iHJ{r)-f{r)Hk) - {H,f{r)-f{r)Hk)Hk+ (6) 
H,iHJ{r)- f{r)Hk) - {HJ{r)-f{r)Hk)Hk)\nklk). 

This is a recurrence relation valid for arbitrary radial potential energy functions, V{r), 
introduced in 

2.2 The Blanchard sum rule 

The second hypervirial takes a particularly simple form when /(r) is a power of the position, 
let us say /(r) = r'*'+^; using this expression for /(r) and restricting ourselves to the Coulomb 
potential, V{r) = —Z/r, we obtain, after a shorter calculation than in |H], the Blanchard 
relation 
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X{Ei - Ekf {nili\r^^^\nklk) = Co{nili\r^\nklk) + Ci{nili\r^ Vfc 4) (7) 

+C2{ni k |r^"^|rafc h)] 

where the hydrogenic energy eigenvalues are Ea = —Z^/2n^, independent of /, and 



Co 

Cl 
C2 



ili-lk)ili + lk + l) 



nt 



nt 



+ A(A + 2) 



-2ZX{\ + 2){2\ + l] 



2 (A + 2) 



(9) 
(10) 



2.3 The Pasternack-Sternheimer selection rule and the non - re- 
lativistic virial theorem 

From this result we can also obtain, as special cases of the relation, the Pasternack - Stern- 
heimer rule jSOl 



{nik \ ^\nklk) = 



which says that the matrix element of the potential vanishes between radial states of 
central potentials when 1) their angular momenta coincide and 2) their energy eigenvalues 
depend only on the principal quantum number. In the completely diagonal case (n^ = Uk, 
li = Ik), we can further obtain the non-relativistic quantum virial theorem|i3J 



{V) = -Z{-)=2{E). 



(12) 



As we exhibit in section 3, we can obtain analogous results using our recurrence relations in 
relativistic quantum mechanics. 



3 The relativistic case 

The method sketched in the previous section can be extended to the relativistic Dirac case. 
To that end, we need to start with the equivalent of the non-relativistic radial Hamiltonian 
of the Schrodinger equation for a central potential V{r) (Eq. 3) in the relativistic case. To 
obtain such expression we begin with the Dirac Hamiltonian Ho for a central potential 

Hd = ca-p + Pc^ + V{r), //^^(r) = ^^(r); (13) 

where we are using again atomic units, a and (3 are the 4x4 Dirac matrices in the Dirac 
representation 
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Here the I's and O's stand, respectively, for 2x2 unit and zero matrices and the a is the 
standard vector composed by the three Pauh matrices a = {a^, Cy, a^)- Since the Hamiltonian 
Hd is invariant under rotations, we look for simultaneous eigenf unctions of Hd, | Jp and J^, 
where J = L + S and 

Hence the solutions of the Dirac equation can be written in the alternative but entirely 
equivalent forms 



(16) 



-Km, 



where Xnm^ and X-Km^, or and 3^j^, are spinor spherical harmonics of opposite parity, 
and K, = —e{j + 1/2) is the eigenvalue of the operator A = /5(1 + E ■ L) which commutes with 
Hd (where H = a® I = diag(cr, a)). Parity is a good quantum number in the problem because 
central potentials are invariant under reflections; parity varies as (—1)' and, according to the 
triangle's rule of addition of momenta, the orbital angular momentum is given by / = j±l/2. 
But, instead of working directly with parity or with k, we prefer the quantum numbers j 
and e, introduced above, which can be shown also to satisfy Z = j + e/2 in all cases. We also 
define /' = j — e/2; accordingly, the spherical spinor depends on / whereas the spherical 
spinor 3^j^, which has the opposite parity, depends on Writing the solutions in the form 
(16) completely solves the angular part of the problem. 



3.1 The relativistic radial Hamiltonian 

To construct the radial Hamiltonian, we use the relation 



{a-r){a-p) = (S-r)(E-p) =r-p + iE-L] (17) 

we then use = [L + (1/2)E]^ = + E • L + 3/4. For the term L ■ E we also need an 
expression for acting on the eigenfunctions (16). Directly from this equation we see that 
when is applied to any central potential state, the big component of the state function 
behaves with the orbital quantum number I = j + e/2, whereas the small one does so with 
the orbital quantum number V = j — e/2; we have then, 

/(/ + l)=j(j + l) + 6(j + ^) + ^, (18) 

for the big component, and 

l'{l' + ^)=3{3 + l)-<J + \) + \, (19) 

for the small one. The action of upon a solution of the form (fTB|) is therefore always of 
the form 

L'=j{j + l)+(3e{j + \) + \, (20) 
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where 13 is the Dirac matrix ()14|) . From this result we obtain the term L ■ E and, substituting 
it into (a ■ p), we finally obtain 



c{a -p) = ar [pr - ipc^-{j + ^)], 

r z 



where 



ar = -a ■ r, pr 



r \ ~^ dr] 



(21) 



(22) 



We are now ready to write the relativistic radial Hamiltonian, and the corresponding 
radial Dirac equation, as Hk%pk{r) = Ekipklr), with given by 



Hi. = car 



Pr - (jk + ^ 

r \ 2 



+ f3c^ + V{r), 



where we introduced the purely radial eigenfunctions 



-1 
-1 



(23) 
(24) 

and the radial Dirac 



in a 2 X 2 representation where, (3 = diag(+l, —1), a, 
equation becomes then jT21 HHI 

+ (Vkir) - Ek) tkC {jk + 1/2) /r - d/dr 

(25) 

_ekc{jk + 1/2) /r + d/dr -c^ + (\4,(r) - E^) J [Gr,,j,e,ir) _ 

We want to remark that though this explicit representation can be used for our problem 
[T3J[T7j, it is not really necessary, since all our results are representation independent. 



3.2 The relativistic hypervirial result 

The relativistic recurrence relation we are after, can be deduced using a similar reasoning as 
the used in section 2 for the non-relativistic case, that is we need a hypervirial result. Let 
us first calculate the non-diagonal matrix element of an arbitrary radial function /(r) 



{E2 - Ei){n2j2e2\f{r)\nijiei) = {n2i 2^2^ 2 f{r) - f{r)Hi\ni jiei) = 

-zc(n2j2e2|ttr (fir) + ^Pfir)j jiCi), 



(26) 



where from now on the labelling in the kets stand for the three quantum numbers n^, jk, 
and ek, we have defined A21 = £2(2^2 + 1) — ei(2ji + 1), and the matrix elements of radial 
functions are calculated as 



{n2j2 e2|/(r)|nu-i e^) = J f{r) {F;{r)F^{r) + G^(r)Gi(r)) dr, 
(ri2j2e2|/3/(r)|nijiei) = j f{r) {F*{r)F^{r) ~ Gl{r)G^{r))dr. 



(27) 
(28) 
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The subscripts stand for the 3 quantum numbers specifying the state. 

We next proceed to calculate a "second order iteration" by substituting /(r) 
H2f{j) — f{r)Hi in the last expression. Let us calculate first and ^-ffi, 



H2i = -c^ {f'/r - f" - f'i - (A2-i/2r/3)/' - iA^j2r)f3fi) + 



(29) 



e2C^ (2j2 + l)/2rP 



tear 



r + 



[V{r) - c^(3] 



and 



A 

dr 



iH, = -(1/r) [c'f - c\A,j2r)Pf\ - [c'f - c\A,j2r)Pf 
-c' {2j, + 1) ^(3 (f - ^(3f) - tear (f + ^f3f] {V{r) + ^c') 



(30) 



and then, we write down the difference of the matrix elements associated with Eqs. and 




{E2-E^f{n2j2e2\f{r)\n,j^ e,) = 
{n2j2e2\[c'{-A,,/2r')Pf -c'f" - c\A,j2r)(3f - c\A,Jr)(3fi+ 
c\At,/2r)Pf' + (A2i/2r)' / + 2^m,/3 (c^ + iA,j2r)Pf)]\n,ji ei). 



(31) 



where we have defined A^^ = e2(2j2 + 1) + ei(2ji + 1). Please notice that here and in what 
follows we are always assuming A21 7^ 0. 

This last expression ()31|) is the direct relativistic equivalent of the generalized second 
hypervirial ^[Cf. Eq. (pTT|) above]. The expression involves the operator d/dr, but here, due 
to the presence of Dirac matrices in the result, we cannot use the trick employed in the 
non relativistic case where we took advantage of the Hamiltonian to simplify the calculation. 
Instead, let us calculate the following second order iteration for non-diagonal matrix elements 



{n2j2 e2\H2i + iHi\ni ji ei) = {El - El){n2]2e2\f{r)\niji ei) = 
(^2j2e2| [-2cV'(r)/r + c2(A2l/2r2)/5/(r) -cV"(0 -2cV'(r)|:+ 
c2(A+ A2l/4r2)/(r) - 2zcaJ/'(r) + ^ (3 f {r)\v {r)] \nai ei). 



(32) 



Due to the presence of Dirac matrices in our results, we also require to calculate non-diagonal 
matrix elements for expressions involving arf{r) and l3f{r), namely 



H2 {-icarf{r)) = 

-c^f/r - c^f - cVf + c\e2/2r) (2j2 + 1) Pf] + ic'arfif - icarV{r)f. 



(33) 



and 
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{-icarf{r))Hi = 

f [c\l/r) (l + r£) + c\ei/2r) (2ji + 1) p] - ic'ar(3f - icarV{r)f- (34) 



adding up these two last expressions, we get 



\n2 



{E2 + Ei){n2j2 62! (-icarfir)) \niji ei) = 
j2 e2\c' (-2//r - f - 2/f ) + c\A^j2r)Pf - 2zcarV{r)f\nai ei). (35) 



(36) 



From the matrix element of H2 {—icarf{r)) — {~icarf{r)) Hi, we can obtain 

-ic{E2 - Ei){n2j2e2\arf{r)\niji ei) = 
(^2^2 621 [-cVW + c2(A+/2r)/?/(r) + 2^cV/?m/(r)] Im ji ei); 

proceeding in a similar way for H2 {[3f{r)) + {Pf{r)) Hi, we get 

iE2 + Ei){n2j2e2\Pfir)\nijiei) = (^2^2 £2! [d/^a^/' - 2car(A2i/2r)/ 

+ 2{c^ + (3V{r))f{r)]\nijiei). 

Equations (j^-ITTjl are the basic equations of our problem. 
3.3 The relativistic recurrence relations 

To proceed from ()37|). we now consider, as in the non-relativistic case, only radial functions 
of the form /(r) = and we insert the explicit expression for the Coulomb potential: 
V{r) = —Z/r. Let us mention though that our results can be generalized to other power 
of potentials, such as the Lennard- Jones potential |^. Substituting /(r) = r'^ in (jSH), it 
follows that 



(37) 



(E| - Ef) (n2j2 (^2\r^\niji ei) 



('^2j2 £2! 



(l/4)c2A+A2l-c2A(A + l) 



„A-2 



+ (cV2)A2l/3r 



A-2 



2c2^^A-l d 
ar 



2tcar A + (1/2) A2I/? r>^''V{r)\niji ei) 



(3^ 



hence, we can eliminate the term containing the derivative operator in this last equation 



using /(r) 



„A-1 



in Eq. to get the result 



(7^2 J2 e2|(c/4)A^iA 



{El - El) (n2j2 ^2\r^\ni3i Ci) 



21' 



,A-2 



+ (c72)A^i/3 (1 - A) r^-'' - icarpA^ir^^^Vir 



+ {E2 + Ei)X {-icar) r^'^\ni ji ei) 



(39) 
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we can, in this last equation, eliminate the term with —icar^2il3r^ ^ by using Eq. (f26|) with 
/(r) = r'^"^, to get 



{El- El){n2j2t2\r^\niiiti) = (722^2 £2! 



> A,, A+ 



+ C 



,2 a: 



il-\)p 



,A-2 



2Z 



icarr^-\l - A) - {E2 - Ei)r^-^ - {E2 + Ei)Azca^r^- Ji ^i) 



Now, from Eq. ()36|) with /( 



r] = r 



A-l 



we get 



(40) 



(^2 - -Ei)(?^2j2e2| -ica^-r^ Vu'iei) = 
(ri2j2e2| - c2(A-l)r^-2^_(c2/2)A+/3r^-2 + 2zcV/5mr^~>iJiei) (41) 



and using /(r) = r'^ in Eq. P7j) to eliminate the term 2ic^arl3mr^ ^ of the above equation, 
we obtain 



{E2 - ^i)(n2j2e2| (-ica^r 
(n2 J2 62! [-c=^ (A - 1) + (cV2)A+ /5r^-2 

+ (cVA) (-ica,) Aslr^-i + c2((4/A)r^ - (4Z/A)/3r^-i 

which can be written as 



^'^^ WiJi ei) - 

- (2/A) {E2 + ^1) c2/5rn 
1^1 jiei); 



(42) 



(^2 - ^1) - c^Asl/Aj (r^2J2e2|(-^ca,r^-l)|r^lJlel) = (^2 J2 e2|[(l - A)r^-2 

+ ^r^ + f A+/3r^-2 - (^)/3r^-i - f (E2 + Ei)/3r^]|nui ei). (43) 
We can also obtain a new relationship for the matrix elements of —icarT'^''^^ using Eq. 

A-l 



(|26|) with /(r) = r^, and substitute the result in Eq. (j4ip to eliminate the term 2ic'ar(3 



mr 



{E2 - Ei){n2j2e2\[-icarr^ ]|?^ijiei) 



{n2j2 e2|[ - c2 (A - 1) + (cV2)A+ /3r^-2 + (4A/A21) (-zm,) r 



-icyA2,)iE2-E,)r^]\naiei). 



(44) 



Rearranging terms, we obtain 



{E2 - El) - ^ 



21 



{n2j2e2\i-icarr^ )|?^ijiei) 



{n2j2e2\c'{l - Xy-' - {E2 - E^y + (f ) A+ /Jr^-^lnui ei). 



(45) 



The relation we are looking for follows from this last result and Eq. (j¥T|l . We use 
succesively r'^~^ and r^"^ from Eq. ()45|1 to eliminate the terms 2(£^2 + Ei)Xicarr^~^ and 
2icarr^"2(l - A) that appear in Eq. to finally get [12] 



Co(^2 j2 e2k^|^i ji ei) 
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(46) 



i=l 



where the numbers Cj, i = 0, . . . 3 are given by 



Co 



Cl 



C2 



C3 



(i?|-i?l^)(ig2-i?l)A21 

(E2-Ei)A2l-4Ac2 ' 
2Z{E2 - Elf 
~(^2-^i)A2l-4(A-l)c2' 
c^A^lA+i _ _ c\Ei + E2)A2-i 



4 " "'(^2-^i)A^i -4Ac2 

2c2Z(A-l)(A-2)A2i 



(E2-i5i)A2-i-4(A-l)c2^ 
and the numbers di, i = 2 and 3, by 



(47) 
(48) 
(49) 
(50) 



do 



d. 



21 



A) 



X{E2 + Ei)A+ 
{E2 - ^1) A21 - 4c2A^ 



Zc\\ 



l)A2iA+ 



(E2-Ei)A2l-4(A-l)c2 



(51) 
(52) 



It seems to be natural that we obtained six matrix coefficients instead of only three as 
in the non-relativistic case. This is to be so, because in the Dirac case we have to deal at 
the same time with the big and the small components in the state function of the system, 
doubling the 'degrees of freedom' we need to know to determine completely a state in the 
relativistic realm. 

Unfortunately it is not easy to avoid the /3-dependency in Eq. (jlBj) . and thus, taken 
on its own, Eq. ()45|) does not allow the computation of < 2 | r'^ | 1 > in terms of the 
< 2 I r^""- I 1 >, a = 1, 2, 3. The situation is not hopeless though, because it is still possible 
to obtain another recurrence relation for non-diagonal matrix elements of (3r^ simply by 
eliminating the term —icarr^~^ between Eqs. (|43|) and (j45j) . In such a way we get 



eo(^2 j2e2|/?r^|ni jiei) = 60(^2 J2 e2k^|^iji ei) + ^2(^2 j2 e2|r^ jici) 
+ei (n2 j2 e2\Pr^~^\ni ji ei) + 62 (n2 j2 £2 l/^r^-^ |ni ji ei) , 

where the numbers 6, and Cj i = 1,2,3 are given by 



(53) 



bo 

Co 

ei 

62 



4 A [(^2 -Eif-ic"^ 
c^(l-A) [{A^,r-4X\ 
2{E2 + Ei)[{E2-Ei)A2i-4Xc''] 
4c^Z[4Xc^ -{E2- El) A21], 
oA^ 



^21 



[(A2-i)^-4A^ 



(54) 
(55) 

(56) 
(57) 

(58) 
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We consider Eqs. fl46|) - (58) as the equivalent of the Blanchard relation for the relativistic 
case. The reason appears to be evident in the next section, but let us mention that these 
equations produce the relativistic virial theorem, as happens to be also the case in the 
non-relativistic Blanchard sum rule (see subsection 2.2). 

To end this subsection, we should say that the relativistic recurrence relations can be, 
juggling with all the relations used (see [TH [T3] for details), put in the uncoupled (and 
possibly the simplest) form 



where 



and 



where 



1^0 



1^1 



1^2 



^3 



1/4 



^5 



J^o{n2j2<^2\r^\nijiei) = ^i^i{n2j2<^2\r^ *|nijiei). 



i=l 



2{E+f{E-fJ 
c^Z{\-l) ' 

-8E+(E 

2AF 



2 J + 6c2 
c2(A-l) 



{X-1)H 
J 



J-E+~ 



Ac\\ - 2) 
S 



E+D 


'G 


2^' 


Z 


~2 


K 



^ ' 2HJ 



c\X-l)L 
2Z 



\E+ 



X ■ 



c^AtiL ( „ Ac 



8Z 



J-E- 



Ac\X - 2) 
S 



8Z{\ - 2)HJ 



K 



^{\-2)L 



Vo{n2j2e2\/3r^\nijiei) = J2^ii^2j2(^2\/3r^ >ijiei). 



1=1 



(59) 



(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

(66) 
(67) 

(68) 



E+D KE+ 2X 



V2 



2XF 
■K 



2& 

D 

XF 



A+ ' 



^21 



2XF 



E+K c2(A-2) 



21 



A, 



21 



4(A-1) 



(69) 
(70) 
(71) 
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(A - 1)ZKL 



(72) 



The symbols defined in order to write the above recursion relations are, M = X{X — 1)E~^, 
D = A^^E~ - Ac^X, F = {E- f - Ac^, ^ = A^^ + A+ , K = ^/A+ , L = AX^ - (A^J^ 
J = (D + 4c2)/A+, G = [J{2M-AtiD)/{X-l)D], and H = D/{D + A^). One more thing 
is worth mentioning, all the relativistic sum rules we have discussed are valid whenever 
the following condition holds u;i + u;2 + |A| + l > where Ua defined by the quantities 
+ \J (ja + 1/2)^ — {ZapY, a = 1, 2 are real numbers, and you must keep in mind that the A's 
are even allowed to be complex [HI 1121 UHl HE] • 



3.4 The diagonal case 

In the calculations of the last section we always assumed that A^^ ^ 0, but some interesting 
results are also obtained when this quantity vanishes. In order to study the diagonal case 
we put ei = €2 and ji = j2', this in turn imply Agi = 0. We remark that in all cases 
A+ ^ A+ ^ 0. 

We start from Eq. (pUj) and put 

{E2- Ei){n2je\f{r)\nije) = (ris j e |(-ica^/'(r) |ni j e ). (73) 

We proceed to calculate the second order iteration by substituting /(r) ^ = H2f{r) — 
f{r)Hi, as was done in the previous section, in the above equation to obtain 

{E2 - E,f{n2j e\f\n,j e) = c^{n2j e\[ - f" + ^f'P + 2zcarPf']\nu e); (74) 
and then substitute /(r) ^+ = H2f{r) + f{r)Hi again in (58) to get 



{El-El){n2je\f{r)\mje) 



(^2 je 



c2(2/'(r)/r + /"(r) + 2/'(r) 



dr ' 



+2icarf' {r)V{r) |nije). 
So, the relevant equations for the diagonal case are now 



(75) 



{E2 + Ei){n2je\{-icarf{r))\nije) 



-(n2j6|[c2(2/(r)/r + /'(r) + 2/(r)|:) 
-2icarV{r)f{r) |nije). 



(76) 



and 



{E2 - Er){n2 3e\{-icarf{r))\nae) = (ns j e |[ - ^f\r) + c2(A+/2r)/5/(r) + 

2tc'ar(3f{r)]\nae). (77) 

We also have, for the matrix elements of f3f{r), 



{E2 + E,){n2je\(3f\n,je) = {n2je\[- tcar^f + 2 (c^ + pV{r)) f]\n,je). (78) 
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We can now obtain a recurrence relation valid in the diagonal case. First, let us use 
/(r) = in Eq. ((711) to get 

{E2 - Eif{n23 e \r^\nij e ) = A(n2 j e |[ - ^{X - l)r^-^ + {c^ /2)^+ Pr^''^+ 

2ic^ar(3r^-^]\niie). (79) 

Evaluating now equation (fTHj) with /(r) = r'^, we obtain 

{E2 + ^1) {n2 3 e \(3r^\ni je) = {n2je\[- icarlSXr^-^ - 2Z (/3/r) + 2cV^] j e ), (80) 
and eliminating the ica^/^Ar^"^ between the last two equations, we finally get 

[(^2 - Ei)2 - 4c^] (n2je|r^|r2ije) = -(cV2)AA+ (^2 j e |/5r^-2|ni j e ) 
-4c^Z {n2j e\l3r^^^\nij e) - 2c^{E2 + Ei) (^2 je | l3r^\nije) 

-c^\{\ - 1) {n2je\r^'^\nije). (81) 

This is the only recurrence relation we get in the diagonal case. 

3.5 The relativistic Pasternack-Sternheimer rule and the relativis- 
tic virial theorem 

The special case when A = in (jHT|) is of particular interest 

[(E2-^l)'-4cVnin. = 

-4Zc2(n2je| (f) Imje) - 2c^iE2 + Ei) {n2j e\P\n,j e). (82) 

What we have obtained in this last expression is the relativistic generalization of the 
Pasternak- Sternheimer rule of non relativistic quantum mechanics jSOl- This rule says that 
the expectation value between hydrogenic states of the 1/r^ potential, vanishes when the 
orbital angular momenta of the states 1 and 2 coincide, that is when h = h- We want to 
remark in this point that in the relativistic case, the expectation value of the 1/r potential 
(which could be regarded as the square root of 1/r^) times j3, does not necessarily vanish 
even when the total angular momenta of the two states coincide: ji = j2- This point 
agrees with the fact that the non-relativistic Pasternack-Sternheimer rule is applicable to 
eigenfunctions of potentials whose energy eigenvalues depend only on the principal quantum 
number — which is not the case for the hydrogen atom in Dirac quantum mechanics [T^ . 

Furthermore, two special cases are immediately deduced from this last expression: 

1) The first case, when rii 7^ n2, is 

2{n2je\^\n,je) = -{E2 + E^){n2J e\P\n,j e). (83) 
r 

2) The other case follows when rii = 77-2 
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c' = - {(3V{r)) + E {P) = z(^^'^+E {(3) , (84) 

which is the relativistic virial theorem [ST]. Also, from the relation < (3 >= [S], we 
can also write an interesting result for the average value of the Coulomb potential and the 
/3 matrix. 

= c^l3V{r)) + c^ = -Zc^/^(^^ + c\ (85) 

4 Analytic results for and /3r^ 

The recurrence relations found above involve expressions that can be burdensome to han- 
dle, excepting perhaps in the diagonal case. Given such situation, before we generalize the 
relations above to the case of twopotentials, we calculate explicit formulas to evaluate the di- 
agonal and the non-diagonal matrix elements of interest in the particular case of the Coulomb 
potential V{r) = —Z/r. The results obtained are based on properties of the hyper geometric 
function and can be deduced directly from the two differential equations that follow directly 
from the Hamiltonian (j23j) . In this section we mostly use ideas from [121 • First of all we must 
say that we are interested in the bound states of the problem, so the quantity k = \J — E'^ 
is positive. We can write the differential equations for the radial part of any central problem 
in terms of the dimensionless variable p = kr. We also define the following quantities 

1 



V^'^W, C = Z/c = ZaF, =e(j + i), 
v^^[c^-E)/{c^ + E), s^^r[^\ (86) 

where ap — 1/137 is the fine structure constant. 

It is proven in that the radial Dirac equation with potential V{r) = —Z/r is com- 
pletely equivalent to 



(+^ + ^ji^(p) =(^-^ + C/p)G(p); (87) 

where F{p) and G{p) are the radial solutions of the big and small components repectively. 
Now, to solve these two equations, we look for solutions of the form 



and 



G{p) = (p) + ^_ (p)] . (89) 
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The solution to these coupled differential equations can be written in terms of the La- 
guerre polynomials of non-integer index |^ ^] 



(90) 



where the Laguerre polynomials C"{p) are related to both the hypergeometric function, 
iFi{—n, a + 1; p), and the Sonine polynomials, T^^\p) jlS], through the relation 



C-n{p) = « + 1; p) = (-l)"r(a + n + 1) T ^(p), (91) 

and a and h are constants. Substitution of these results in Eq. (jHTj) gives the condition 



a{Tj + s- + n)+ b{n + 2s) = 0, 
b{Tj — s + Qu^^ — n) — an = 0. 



(92) 



When we solve the above equations we obtain a relationship between n and v. FromEq. 
(j86|l we see that this in turn gives us a an expression for the energy, whcih coincides with 
the correct value provided we define the principal quantum number N = j + l/2 + n, where 



n = 0, 1, 2, 



E = mc^ 



1 + 



\ (iV-j - 1/2 + ^(j + 1/2)2-^2 



To proceed further, we take 

h = -a{Tj + s + n- Cv"^)/{n + 2s), 
and write the result in a symmetrized form: 



-1/2 



(93) 



(94) 



F(p) = Vmc2 + E Cp'e-P [u Cl'{2p) + v Cl'_^{2p) 
G{p) = -y/mc^ - ECp'e-P \uCl'{2p) - v Clt^{2p) 



where 



u 



{vj + s + n- C,v 



-lNl/2 



(95) 



(96) 



The C in Eq. (|95|) is a normalization constant that can be obtained from 

/■°° V(n n ^^ 

JO n\ 
To obtain C, we use relations (jHSI) to get {Tj + s + n — Cu^^)'^ = {n + s — Tj — (h'~^)/n{n + 2s); 
we need also (n + s) = (E/y/m'^c^ — E"^, which is obtained from the expression for the energy 
eigenvalues of the Dirac hydrogen atom: after some work we obtain 
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The calculation of diagonal, arbitrary power matrix elements of the form < r'^ > and 
< /5r^ >, are obtained from 

InmW^ e--x"+^£^(x)C(x)rfx. (99) 
Jo 

This expression converges for Re{a+\+l) > 0, and is zero if A is an integer such that m—n > 
A > 0, where without loss of generality, we assume that m > n. From Rodrigues formula 
and {dJ^ / dx'^)x^^^ = (—!)'"[— A; — X\mX^'^^~"^ , where [n], n an integer, is a Pochhammer 
symbol, we find, after a m-times partial integration, 

C.(A) = 1- + „ + i)r(„ + t + A + l)[-k - A|,„ ^ 

ml k\{n-k)\T{a + k + l) ^ ' 

We use now the identity \—k — A]m = \—k — A]fc[— A]m-A;, change the order of summation 
k ^ n — k and use the identities 



A]m-n+fc — [— A]m_„[— A + m — 

T{n + a + l) = {-l)''r{a + n- k + l)[-a-n]k, 

T{a + X + n + 1) = (-l)^r(a + A + n - A; + l)[-a - A - njfc, 

i.-,,,-A]„_. = (-i)" ''r(i:y'i_,^„i, . (101) 



to deduce that 



r(a + A + n + l)r(A + n + 1) 



m!n!r(A + 1) 

X3-^2(— — n, —A + m — n, — n; —A — n, —a — X — n; 1). (102) 



Now, we are going to consider two cases for the matrix elements (n2 j2 ^2!^ I^i Ji ^i) and 
(77,2^2 ^2\P^^\ni ji ei); in the first one ki = k2, where we need to evaluate 

KTW^ x'^^'^^\-^C^t'\x)C^^'^\x)dr, (103) 
Jo 

and the second one, when ki k2, where we need 

Km'W= r''+'''+^e-^'''+^'>C^^''\2kir)C^^''\2k2r)dr. (104) 
Jo 

In the first case, we see that integral p03|) is convergent if Re{si + S2 + A + 1) > 0, and 
vanishes when si — S2 + A is an integer such that m — n > si — S2 + X > 0. With a similar 
reasoning as in the diagonal case, we get 

A-r=(A) = H + + A + „ + l)r(A + + „ + 



Im—n 



ml n\ r(A + Si — §2 + 1) 
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X 3-F2(— 2si — n, w + m — n, — n; w — n, —A — Si — S2 — n;l). 



(105) 



where w = —A + S2 — si. In the second case we see that the integral converges for Re{si + 
S2 + A + 1) > 0, and is not zero provided ki k2- A calculation by parts shows that 

r.s,s,(,. _^^ {-iy{k2 - fci)— (fci + A;2)(— ^-^-^-^) 

T{n + 2si + l)[s2-si-X-j]i , \ ■ , IN nnfiN 

X ^ ■ ^ n ^ + S2 + si + A - z + 1 , 106 

r(2si + J + 1) 

where ki ^ k2. This expression can be simplified a little bit by eliminating the double sum 
in terms of a single one using the identities. 

[p]m-^ = i-ir-'j ^l^rf^^ TTV (107) 

[—p — m + IJj 1 {—p — m + 1) 
T{-p-m + i + l) = [-p-m+l]iT{-p-m + l), (108) 

(m-z)! = -j^, (109) 
T{p + 1) = [p-m + l]^T{p-m+l). (110) 

This expression is validad provided p, m and i are integers. After some algebra, we finally 
get 



^nm K"^) m! (fe2+fci)''i+'2+^+i ^j=0 j! r(2si+j+l) 

X2Fi(-m, Si + S2 + A + 1; Si - S2 + A + j - m + 1; ^^^^ 



fe2+fcl 



;iii) 



We are now in the position to calculate the values of < r'^ > and < (3r^ >, using the 
expression for J"^(A), with a = 2s. After some algebra we obtain from Eq. and 
our previous definitions the following expressions 



and 



< pr^ >-- 



mc^\C\'' 
(2A;)^+i2^-i 



{2kY+^2'-^ 

The numbers u and v are the constants defined in Eq. 



(112) 



(113) 



5 Generalized recurrence relations for two potentials 

In this section we exhibit that the hypervirial techniques suffice to obtain relations between 
relativistic matrix elements of an arbitrary radial function between eigenstates corresponding 



18 



to two different potentials. We obtain four recursion relations between the matrix elements of 
arbitrary radial functions, /(r) and /5/(r), with the matrix elements of their first and second 
derivatives, taken between eigenstates of different potential functions, Vi(r) and V2(r) — both 
behaving as the fourth component of a 4-vector. We call these relations between matrix 
elements generalized recurrence relations because they led to relations relating succesive 
powers of the radial coordinate when the arbitrary function /(r) is chosen to be where A 
is a constant number. The recursions between matrix elements we deal with in this section 
are calculated under the further assumption that the two potentials, Vi(r) and V2(r), have 
a common minimum, that is, that they are referred to the same origin. The relations lead 
to explicit recursions when specific forms of the radial functions are used and particular 
potentials are substituted. The recursions can be useful for studying radiative transitions in 
Rydberg atoms, in analysing atomic photorecombination and photoionization processes, for 
example, any transition to an autoionizing state studied in the central field approximation 
where the electron motion is solution of the Dirac equation with an effective central potential 
created by a or a — x electron ion core; or in any other atomic processes involving highly 
excited electrons which need to be studied using multichannel spectroscopy or quantum 
deffect theory [2SI 113 HE! US] • It can be also useful for calculating relativistic corrections 
to ionic oscillator strengths, or in analysing impact ionization or vibrational transitions in 
molecules — albeit in the last two cases in a rather crude manner jSOl |23 EH IM| • 

Let us mention that analogous relations for matrix elements between eigenstates of two 
potentials were obtained some years ago in nonrelativistic atomic physics also with the help 
of hypervirial results [HI EH] ■ 

Let us consider two radial Dirac Hamiltonians with two possibly different radial potentials 
(each behaving as the temporal component of a fourvector) Vi{r) and V2(r). We further 
assume that these potentials have the same equilibrium position which, furthermore, is 
coincident with the origin of coordinates. That is, the recurrence relations correspond to the 
so-called unshifted case |9^. The main difficulty for not dealing with the general (so called, 
shifted) case is the angular couplings introduced by the relative displacement of one of the 
potentials respect the other in the otherwise purely radial interactions. 

These two Dirac Hamiltonians can be written as 

Hi = car[pr - ^/3ei(ji + l/2)/r] + Mi/Jc^ + V,{r), (114) 

and 

H2 = car[pr - iPe2ij2 + l/2)/r] + + V2ir), (115) 

where we are assuming Mi 7^ M2 as it is convenient for our calculations. The eigenstates 
of these radial Hamiltonians correspond to a definite value of the total angular momentum 
J = L + S and of the quantum number e. The Dirac equations are Hkipkir) = E^tpklr) where 
the energy eigenvalues E^^j^g^ = and the corresponding eigenfunctions i^n^j^Ski.'^') = i^ki^) 
are assumed known. 
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5.1 The first generalized recurrence relation for the case of two 
potentials 

Taking the difference between the radial Hamiltonians Hi and H2 in ()114|1 and we 
obtain 

Hi = H2 + icar(3^ - c^(3M- - (V^2(r) - Ki(r)) . (116) 
2r 

where = M2 ± Mi , and A"*^ = 62 (2j2 + 1) ± Ci (2ji + 1) and, in general, if X is any symbol 
we define = X2 ± Xi. We can directly evaluate the commutator 

[Hi,fir)] = -zcar^ (117) 
ar 

where /(r) is an arbitrary radial function and [H, /(r)] stands for the commutator between 
H and /(r). We can calculate this commutator again, but now using equation ()116p . to get 
the alternative form 

[Hu f{r)] = H^fir) - f{r)H, + (^zca,/5^ - c^/JM" - V'^ f{r). (118) 
It is now simple to obtain, from equations ()117|1 and ()118|1 . the relation 



{E2 - i?i)(2|/|l) = (2| {c'[3M- + V) f\l) - ic{2\ar [f'' + /?^/^ |1); 



(119) 



where we have additionally taken matrix elements between the eigenstates (1| = (nijicil 
and \n2j2^2) = |2), and we have defined 

= V2{r)±Vi{r). (120) 

Equation pi9p leads to recursions between relativistic matrix elements of radial functions 
between hydrogenic states ( [SI CHI EH] ) , and generalizes a nonrelativistic expression useful 
for similar purposes |6^. On the other hand, it is an exact relation for the calculation of 
any matrix elements of /(r) between eigenstates of two potentials in relativistic quantum 
mechanics. 

Taking the potentials as equal, i. e. Vi{r) = V2{r), we recover a relation which has 
been useful for obtaining recurrence relations between atomic matrix elements in relativistic 
quantum mechanics [T3|[Tn]. Albeit exact, equation (|119|) is not entirely convenient due to 
the presence of the operator ar(3. To get rid of this term, we found it convenient to work 
directly in terms of operators and not in terms of the matrix elements themselves. The 
matrix elements will be evaluated at the end of the operator calculations. 

Let us first establish that 

H2f - fHi = {c^l3M- + V')f- tear (j' + f3f^^ , (121) 
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notice that equation pi9|) above can be obtained from p21|) just by taking matrix elements. 
The following result is also easily established 

H2f + fH, = {c'(3M+ + r+) / - tear (^2/^ + /' + ^ + /?/^) • (122) 
Then, it can be seen that 



- ic{H2arf + arfHi] 



ica 



-c' hf^ + f' + ^-pf ^ 
\ ar r 2r 



(123) 



and that 



H2fV- - fV"H^ = [c^l3M' + V~) f 



(124) 



It is also readily apparent that 



/ / 

— ic\H2arP — h arP—Hi] = —icar 
r r 



((3V+ - c^M 

Let us define tp{r) = H2f{r) + f{r)Hi^ and evaluate 



p 



7' f\ A+/ 



2r r 



(125) 



- = c'p^f + (^) / + {^'f^M' +v-yf 



A 



— zcar 



V-f + ^f + c'M-^f + V-^(3f 
ar 2r 2r 



(126) 



In this way, working with all the previous equations, in the way we have exhibited in the 
previous sections, we can get 



H2^ - ^H, = -c' (^r - f3f^^ + 2c'^Pf + {M-fc'f 



- c'M-V+Pf + V-[2 {H2f - fH,) -V']-c 



,A2 
2r 



Pf 



A+ 
'2r 



■f 



(127) 
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Evaluating the matrix elements between the Dirac eigenstates (2| and |1) and rearranging, 
we finally obtain the relation 

ao(2|/|l) + a2(2|^|l) - 2E-{2\V-f\l) + (2|(l^-) + c'{2\f"\l) = 
bo{2\(3f\l) + bi{2\(3^\l) - c'M-{2\V+(3f\l) 
+ c'M+{2\V-i3f\l) + 64(21/3^11), (128) 



where 



ao= (E-f -(c^M-)'' 
-,2 



--A-A+ 
4 



02 

bo = (m-E^ - M^E-) 

&4 = y (A+ - A-) (129) 



This is the first relation between matrix elements of an arbitrary radial function /(r) be- 
tween eigenstates of two different potentials as a function of the eigenenergies in relativistic 
quantum mechanics. 

5.2 More generalized recurrence relations for the case of two po- 
tentials 

Given that the radial eigenstates have two components in relativistic quantum mechanics, it 
should be clear that we need more relations. We can obtain at least three more, following a 
path similar to the one outlined above and in previous papers to get 

co(2|/|l) + a2(2|^|l) - E+{2\V-f\l) - E' {2\V+ f\l) + {2\V+V- f\l) 

-c^(2|^|l) + c\2\f\l) = |(2|/3^|1) + 64(21/3^11), (130) 

where the only newly defined coefficient is 

Co = E+E' - c^M^M^. (131) 

We also get 

eo(2|/|l) = go{2\(3f\l) - (2| (y+ - V') /3/|l), (132) 
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where 



eo 
9q 



(m+ - M~ 



(133) 



This is a very simple equation that, besides, allows writing the matrix elements of / in terms 
of those of f3f. To take advantage of this fact, substitute equation ()132j) into p28|) to obtain 
the new relation 



Ao{2\f3f\l) + Ai(2|/?4|1) + A,{2\V-(3f\l) + ^3(2! (vY f3f\l) 



+ (2| {v-fpfll) + A,{2\V+(3f\l) + 2E-{2\V-V''Pf\l) 



(2| (V-) V+/?/|l) + a,{2\ (r+ - V-) - c'{2\ (V+ - V-)"pf\l) 



A,{2\pl^\l) + 2c'{2\ (y+ - V-)'pf\l) - c'go{2\Pf"\l) 
-c'{2\lv+-V-)pf"\l). 



(134) 



where the newly defined coefficients are 



^0 
Ai 

A2 
A3 
A5 



E-y iE+ -E- 



.2 TP- 



c'E 



,2 / , \ 2" 

M" + ( M 



^M+M- (e+ + E- 



-— [E^ - E-) A+A- - c^A- (M+ - M- 

-2E- (e^ - E-) + {^E-y - (c^M-y - c^M+ (m+ - M" 
E+ - 3E- 



A^ = %- (m+ - M-) (a+ - A- 



(135) 



Equation ()134p is the fourth recurrence relation for the calculation of relativistic /(r) matrix 
elements between states of different potentials in terms of their energy eigenvalues. Notice 
that, at difference of the previous relations [equations (jl28p . (jl3(Jj) and (jl32|) ]. equation (jl34j) 
relates among themselves matrix elements of /9/ and its derivatives times a certain function 
of r. 

In summary, in this section we have obtained generalized recurrence relations for the 
calculation of matrix elements of a radial function between states of two different radial 
potentials sharing a common origin. The obtained sum rules are given in the most general 
case of an arbitrary function taken between any non-necessarily diagonal radial eigenstates 
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of the two radial potentials. Such relations have, as particular cases, sum rules between one- 
potential integrals or, in other particular cases, between overlap and one centre integrals in 
Dirac relativistic quantum mechanics. We expect the obtained formulas, together with the 
previous relations we have obtained [121 El US] , to be useful in atomic or molecular physics 
calculations as they may simplify calculation of matrix elements in the range of applicability 
of Dirac's relativistic quantum mechanics |5H I19j. For most uses of the relations we first 
have to set Mi = M2, i. e. M~ = and M"*" = 2 — assuming the particles are electrons — 
since the use of unequal masses is basically a recourse of our calculational method. 

From a more practical angle, there is little that can be done for the analytical evaluation 
of these integrals beyond the Coulomb and the few similarly exactly solvable potentials. 
However, there are relativistic theoretical techniques [2^1 and numerical methods that, after 
being adapted to relativistic conditions, can provide the crucial "seed" results needed for 
the systematic use of the recurrence relations obtained here |^. Our results can be also 
useful in the so-called perturbation theory of relativistic corrections, in relativistic quantum 
deffect calculations, and for the relativistic extension of the calculations of exchange integrals 
using Slater orbitals or Coulomb-Dirac wave functions [231 IMl IHll liHl 111] • On the other 
hand, our results may be also of some interest in nuclear studies since the 3D Woods-Saxon 
potential, used in the Dirac equation for describing the interaction of a nucleon with a heavy 
nucleus, has been recently explicitly solved and its eigenfuctions can be expressed in terms 
of hypergeometric functions |60j . 

Acknowledgements 

This work has been partially supported by PAPIIT-UNAM (grant 108302). We acknowledge 
with thanks the comments of C. Cisneros and I. Alvarez. We also want to thank the friendly 
support of Dua and Danna Castillo- Vazquez, and the cheerful enthusiasm of G. R. Inti, F. 
A. Maya, P. M. Schwartz (Otello), P. O. M. Yoh, H. Kranken, A. S. Ubo, G. Sieriy, M. 
Chiornaya, P. A. Koshka, G. D. Abdul, F. C. Sadi, and D. Gorbe. 

References 

[1] H. A. Kramers, Quantum mechanics, North Holland, Amsterdam, (1957). 

[2] P. W. Atkins, Molecular Quantum Mechanics, Oxford, Clarendon, (1970). 

[3] O. L. de Lange and R. E. Raab, Operator Methods in Quantum Mechanics, (Clarendon, 
Oxford, 1991). 

[4] Kobus, J. Karkwowski, and W. Jaskolski, J. Phys. A: Math. Gen. 20, 3347, (1987). 

[5] B. Moreno, A. Lopez-Pineiro and R. H. Tipping, J. Phys. A: Math. Gen., 24, 385, 
(1991). 

[6] H. N. Nunez- Yepez, J. Lopez-Bonilla, and A. L. Salas-Brito, J. Phys. B: At. Mol. Opt., 
28, L525, (1995). 



24 



[7] H. N. Nunez- Yepez, J. Lopez-Bonilla, D. Navarrete, and A. L. Salas-Brito, Int. J. 
Quantum Chem. 62, 177, (1997). 

[8] P. Blanchard, J. Phys. B: At. Mol. Phys., 7, 993, (1974). 

[9] Fernandez F. M. and Castro E. A. 1987 Hypervirial Theorems (Berlin: Springer), (1987). 

[10] Fernandez F. M. and Castro E. A. 1996 Algebraic Methods in Quantum Chemistry and 
Physics, (Boca Raton: Chemical Rubber Company), (1996). 

[11] J. J. Peiia, P. Portillo, G. Ovando, and V. Gaftoi, Int. J. of Quantum Chem., 65, 205, 
(1997). 

[12] Martmez-y-Romero R. P., Nunez- Yepez H. N., and Salas-Brito H. J. Phys. B: At. Mol. 
Opt. Phys., 33, L367, (2000). 

[13] Martmez-y-Romero R. P., Nunez- Yepez H. N., and Salas-Brito A. L. J. Phys. B: At. 
Mol. Opt. Phys., 34, 1261, (2001). 

[14] Martmez-y-Romero R. P., Niifiez- Yepez H. N., and Salas-Brito A. L. J. Phys. B: At. 
Mol. Opt. Phys., 35, L71, (2002). 

[15] Ilarraza-Lomeli A.C., Valdes-Martinez M.N., Salas-Brito A.L., Martinez-y-Romero R.P., 
Nunez- Yepez H.N. Int. J. Quantum Chem. 90, 195-201, (2002). 

[16] R. P. Martinez- y- Romero, A. L. Salas-Brito and J. Saldana-Vega, J. Math. Phys., 40, 
2324, (1999). 

[17] Constantinescu and E. Magyari, Problems in Quantum Mechanics, (Oxford, Pergamon, 

1971) . 

[18] I. P. Grant in G. W. F. Drake Editor, Atomic, Molecular and Optical Physics Handbook, 
(American Institute of Physics, Woodbury, 1996) Ch. 32. 

[19] R. E. Moss, Advanced Molecular Quantum Mechanics, (London, Chapman and Hall, 

1972) . 

[20] Bessis N, Bessis G, and Roux D., Phys. Rev. A , 32, 2044, (1985). 
[21] Brack M., Phys. Rev. D , 27, (1950). 

[22] Nana Engo S G, Kwato Njock M G, Owono Owono L C, Lagmago Kamta G and 
Motapon O., Phys. Rev. A, 56, 2624, (1997). 

[23] Owono Owono, L. C., Kwato Njock M. G. and Oumaro B., Phys. Rev. A, 66, 052503, 
(2002). 

[24] Shabaev V. M. J. Phys. B: At. Mol. Opt. Phys. 24, 4479, (1991). 

[25] Wong M. K. F. and Yeh H-Y. Phys. Rev. A 27 2300, (1983), Wong M K F and Yeh 
H-Y., Phys. Rev. A, 27, 2305, (1983). 

25 



[26] Dobrovolska .1 V. and Tutik R. S., Phys. Lett. A 260, 10, (1999). 

[27] Quiney H. M., Skaane H. and Grant 1. P. J. Phys. B: At. Mol. Opt. Phys., 30, L829, 
(1997). 

[28] Elboudali F. and Joulakian B. , J. Phys. B: At. Mol. Opt. Phys., 34, 4877, (2001). 

[29] Ilarraza-Lomcli A. C, Valdes-Martmez M. N., Salas-Brito A. L., Martmez-y- Romero R. 
R, and Nunez- Yepez H. N. , Int. J. Quantum Chem., 90, 195, (2002). 

[30] S. Pasternack, R. M. Sternheimer, J. Math. Phys., 3, 1280, (1962). 

[31] Y. S. Kim, Phys. Rev. 154, 17, (1967). 

[32] West J. B. J. Phys. B: At. Mol. Opt. Phys. 34, R45, (2001). 

[33] Avram N. M. and Draganescu Gh. E. Int. J. Quantum Chem. 65 655, (1997). 

[34] Charro E., Lopez-Ferrero S., Martin I. J. Phys. B: At. Mol. Opt. Phys. 34, 4243, (2001). 

[35] Charro E. and Martin I. J. Phys. B: At. Mol. Opt. Phys., 35, 3227, (2002). 

[36] Nagy L and Vegh L. Phys. Rev. A, 46, 284, (1992). 

[37] Schippers S., Miiller A., Ricz S., Bannister M. E., Dunn G. H., Slachter A. S., Hinojosa 
G., Cisneros C., Aguilar A., Covington A. M., Gharaibeh M. F., and Phaneauf R. F. 
Phys. Rev. A, 67, 032702, (2003). 

[38] Schippers S., Miiller A., Ricz S., Bannister M. E., Dunn G. H., Bosek J., Slachter A. S., 
Hinojosa G., Cisneros C, Aguilar A., Covington A. M., Gharaibeh M. F., and Phaneauf 
R. F. Phys. Rev. Lett., 89 193002, (2002). 

[39] Miiller A., Phaneuf R. A., Aguilar A., Gharaibeh M. F., Schlachter A. S., Alvarez I., 
Cisneros C, Hinojosa G., and McLaughlin B. M. J. Phys. B: At. Mol. Opt. Phys., 35, 
L137, (2002). 

[40] Covington A. M., Aguilar A., Covington I. R., Gharaibeh M. F., Shirley C. A., Phaneauf 
R. F., I Alvarez, Cisneros C, Hinojosa G., Bozek J. D., Domingucz I.. Sant'Anna, 
Schlachter A. S., Berrah N., Nahar S. N., McLaughlin B. M. Phys. Rev. Lett. 87, 
243002, (2001). 

[41] Aguilar A., West J. B., Phaneuf R. A., Brooks R. L., Folkmann F., Kjeldsen H., Bozek 
J. D., Schlachter A. S., and Cisneros C. Phys. Rev. A 67, 012701, (2003). 

[42] U. Fano and L. Fano, Physics of atoms and molecules. University of Chicago, Chicago, 
(1972). 

[43] R. P. Martinez-y- Romero, A. L. Salas-Brito and J. Saldaiia-Vega, J. Phys. A: Math. 
Gen. 31 L157 (1998). 

[44] L. Davies Jr. Phys. Rev., 56, 186, (1939). 



26 



[45] W. Magnus and F. Oberhettinger, Formulas and Theorems for the Special Functions of 
Mathematical Physics, Chelsea, New York, (1949). 

[46] J. D. Bjorken and S. D. Drell, Relativistic Quantum Mechanics, New York, (1964). 

[47] Mannervik S., DeWitt D., Engstrom L., Lindberg J., Lindroth E., Schuch R., and Zong 
W. Phys. Rev. Lett, 81, 313, (1998). 

[48] Xiao-Min Ton, Lei Liu, Jia-Ming Li Phys. Rev. A, 49, 4641, (1994). 

[49] Aymar M, Greene C. H., Luc-Koenig E. Rev. Mod. Phys., 68, 1015, (1996). 

[50] Schippers S, Miiller A, Ricz S, Bannister M E, Dunn G H, Bosek J, Slachter A S, 
Hinojosa G, Cisneros C, Aguilar A, Covington A M, Gharaibeh M. P., and Phaneauf 
R .F. Phys. Rev. Lett, 89, 193002, (2002). 

[51] Semenov S K and Ciierepkov N A 2003 J. Phys. B: At Mol. Opt Phys., 36, 1409, 
(2003). 

[52] Week P., Fojon O. A., Hanssen J., Joulakian B., and Rivarola R. D. Phys. Rev. A, 63, 
042709, (2001). 

[53] Bauche-Arnould C, Bauchc J., Luc-Koenig E., Wyart J.-F., More R. M., Chenais- 
Popovics C, Gauthier J.-C, Geindre J. -P., and Tragin N., Phys. Rev. A, 39, 1053, 
(1989). 

[54] Bang J. M. and Hansteen J. M., J. Phys. B: At Mol. Opt Phys. 35, 3979, (2002). 
[55] Morales J., Phys. Rev. A, 36 4101 (1987). 

[56] Ley-Koo E., Jauregui R., Gongora A., and Bunge C. F. Phys. Rev. A, 47, 1761, (1993). 
[57] Chen Z., Bessis D. and Msezane A. Z. Phys. Rev. A 47, 4756, (1993). 
[58] Kuang J. and Lin C. D. J. Phys. B: At Mol. Opt Phys., 29, L889, (1996). 

[59] Rutkowski A. Phys. Rev. A, 53, 145, (1996). 

[60] Jian-You G., Xian Cheng P., and Fu-Xin X., Phys. Rev. A, 66, 062105, (2002). 



27 



